IMMEDIATE BUCKLING MODEL, HYSTERESIS MODEL, AND CLOTH SIMULATION 
METHOD BASED ON THE INVENTED MODELS, AND COMPUTER-READABLE 
MEDIA STORING A PROGRAM WHICH EXECUTES THE INVENTED SIMULATION 
METHOD 

5 FIELD OF THE INVENTION 

The present invention relates to a method for simulating the movement of 
cloth and computer-readable media that stores a program which executes the 
method, and more particularly to an immediate buckling model, hysteresis 
model, and cloth simulation method based on the invented models, and 

10 computer-readable media storing the program which executes the invented 
simulation method. Deformation by compressive force applied at two extremities 
of the deformation unit which models cloth is simulated by an immediate 
buckling model that is immediately bent without contraction. Deformation by 
stretching force is simulated by a spring model, while cloth hysteresis 

15 characteristics are modeled by using new hysteresis models called "Kwang-Jin 
Choi and Hyeong-Seok Ko Discrete Hysteresis Model " (hereinafter referred to as 
Choi-Ko Discrete Hysteresis Model ) and " Kwang-Jin Choi and Hyeong-Seok Ko 
Continuous Hysteresis Model " ( hereinafter referred to as Choi-Ko Continuous 
Hysteresis Model ), thereby improving the simulation of the movement of cloth 

20 and solving the buckling instability problem. 

BACKGROUND OF THE INVENTION 

In general, a thin structure is characterized as having weak resistance to 
bending force and strong resistance to stretching force, such that movement of a 
thin structure is solved by a stiff equation (a type of differential equation). When 
25 compressive force is applied to a thin structure, contraction occurs. As the 
compressive force is increased to pass a critical point, the thin structure is 
vertically bent, also known as buckling. Buckling is a deformation which occurs 
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abruptly so that it is a very unstable reaction. Therefore, the simulation of that 
kind of buckling results in a divergence problem in differential equations for 
simulating the movement of cloth. 

The buckling instability refers to the contrast of the decrease in 
compressive force and increase in bending deformation. This structural 
instability makes the system matrix extremely ill-conditioned or indefinite, and in 
a case where the time period of the simulation steps is increased, the system 
matrix becomes divergent. In conventional simulations of the movement of cloth, 
a model in which the instability of buckling exists is used such that it is very 
difficult to simulate a phenomenon of wrinkles (buckling) forming on the surface 
of the cloth. 

The buckling instability problem arises not from the stiff equation itself 
but from the structural instability of Jthe cloth. Therefore, simply employing an 
implicit method cannot solve the problem. Resolution of the buckling instability 
15 problem was conventionally sought by adding a damping term to a system 
matrix. However, although the addition of the damping term can stabilize the 
system, it obstructs the naturalness of the movement of the cloth. In other 
words, the damping force prevents formation of wrinkles on the cloth's surface, 
and prevents wrinkles from disappearing. Hence, a simple addition of a damping 
20 term is not desirable to simulate the movement of cloth. 

Furthermore, a continuous body model used in the past brought forth an 
undesirable result due to the following reasons. A coarse discretization should be 
allowed to guarantee a reasonable performance in computer graphic applications. 
However, a continuous body model requires a very fine mesh in order to simulate 
25 large deformation of cloth. Therefore, a reasonable processing speed cannot be 
obtained. Another drawback is that the continuous body model cannot properly 
deal with the divergence derived from buckling, thereby requiring additional 
calculations. [ Eischen et al. 1996 ] uses a non-linear shell model for cloth 
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simulation and performs a finite element analysis to simulate buckling. However, 
great care and cautious measures have to be undertaken such as arc length 
control to prevent divergence due to the non-linearity of the load-deflection curve 
or the singularity of the instantaneous stiffness matrix caused by buckling. 

5 In [ Baraff and Witkin 1998 ], a system of connected triangles was 

proposed as a model. The in-plane deformation energy of each triangle is 
derived from the continuous mechanism. The bending deformation is based on 
the angle between adjacent triangles. Therefore, the bending deformation and 
the in-plane deformation are modelled independently, which allows a large 
10 bending deformation between the triangles regardless of the in-plane rigidity of 
each triangle, thereby resolving the drawback of the continuous body model. 
However, the buckling instability still remains in this model because each triangle 
is modeled as an almost incompressible material and the bending rigidity 
between triangles is very weak. 

15 Still furthermore, in order to simulate the hysteresis characteristics of 

the movement of cloth, the results. of measuring the hysteresis characteristics of 
the cloth have been directly used. However, a hysteretic response curve obtained 
by the measurement can be applied only when a curve change during the 
measurement and a curve change during the simulation are the same, such that 

20 it is not generally appropriate for dynamic simulation. In other words, a 
hysteretic response curve is a function of the entire curvature history, such that 
it is not reasonable to simulate by measuring the curve instead of by 
constructing a physical model for the hysteresis characteristics. 

25 SUMMARY OF THE INVENTION 

It is an object of the present invention to provide a method for simulating 
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the movement of cloth realistically without the requirement of very fine meshing 
and short simulation intervals by using a model without the buckling instability 

It is another object of the present invention to provide a method for 
simulating the movement of cloth and computer-readable media storing a 
5 program which executes the method to accurately express the characteristics of 
the movement of cloth and to achieve a fast simulating speed by solving the 
buckling instability. 

BRIEF DESCRIPTION OF THE DRAWINGS 

10 FIG. 1 is a drawing for illustrating a buckling phenomenon of a rigid 

material; 

FIG. 2 is a drawing approximating a bending deformation of cloth 
according to an immediate buckling model; 

FIG. 3 is another drawing approximating the bending deformation of 
15 cloth according to an immediate buckling model; 

FIG. 4 is an example of cloth deformation in an equilibrium state; 

FIG. 5 is a graph for illustrating the distance between particles positioned 
at two extremities of a deformation unit in relation to the compressive force in an 
immediate buckling model; 
20 FIG. 6 is an approximate graph for illustrating the distance between 

particles positioned at two extremities of deformation unit relative to the 
compressive force at an immediate buckling model; 

FIG. 7 is a graph for illustrating an actual hysteresis response curve 
measured on a pure cotton; 
25 FIG. 8 is a graph for illustrating the construction of a slip unit and a 

response thereof; 

FIG. 9 is a schematic drawing for illustrating a construction that 
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combines two springs and a slip unit, and a response thereof; 

FIG. 10 is a schematic drawing of an embodiment of a construction that 
combines slip units and springs for modeling hysteresis characteristics of cloth 
according to the present invention, and a response thereof; and 
5 FIG. 11 is a schematic drawing of an embodiment of a construction 

combining infinite numbers of springs and. slip units for modeling hysteresis 
characteristics of cloth according to the present invention, and a response 
thereof. 

10 DETAILED DESCRIPTION OF THE INVENTION 

In order to accomplish the aforementioned objects, the present invention 
proposes to use the immediate buckling model, the " Choi-Ko Discrete Hysteresis 
Model " and the " Choi-Ko Continuous Hysteresis Model The immediate 
buckling model comprises deformation units which model cloth, wherein the 

15 deformation unit is bent immediately without contraction when compressive 
force is applied to two extremities of the deformation unit. 

The hysteresis phenomenon of cloth is modeled by using finite numbers 
of spring-slip units in the "Choi-Ko Discrete Hysteresis Model" which is used for 
simulating the hysteresis characteristics between the curvature and moment. The 
20 hysteresis phenomenon of cloth is also modeled by using infinite numbers of 
spring-slip units in the "Choi-Ko Continuous Hysteresis Model" which is used for 
simulating the hysteresis phenomenon between curvature and moment more 
realistically. 

The spring model simulates the deformation of cloth when stretching 
25 force is applied to both extremities of the deformation unit (when stretching 
force is applied to the cloth). Therefore, the movement of cloth in relation to 
compression and stretching that react on the cloth will be analyzed by using 
models respectively appropriate thereto. 
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(ANALYSIS OF BENDING DEFORMATION BY IMMEDIATE BUCKLING MODEL) 



In the present invention, when compressive force is applied to cloth, an 
immediate buckling model (which does not resist to compressive force and does 
not shrink but bends immediately right after the application of the force) is used 
5 to simulate the bending deformation of the cloth. The bending deformation 
energy accumulated on the cloth can be calculated by estimating the shape of 
the deformed cloth. 

Referring to FIG. 2, the process of obtaining deformation energy in 
relation to the bending deformation of cloth caused by the compressive force will 
10 be described below. According to the immediate buckling model illustrated in 
FIG. 2, the region of a predetermined deformation unit of cloth is modeled as a 
straight beam having no surface area and also a plate having a surface area, 
which will be described later. 



15 applied, the predetermined deformation unit is a straight beam of length LI. 
After buckling, the beam eventually reaches a stable equilibrium state (bent 
state). For example, in the present invention, the equilibrium shape is 
approximated as a circular arc with arc length LI. In other words, there is no 
contractive deformation in this immediate buckling model. That is, only bent 

20 deformation occurs. The bending deformation energy can be calculated by 
estimating the shape of the deformation unit using the following equation: 



Where, M is the bending moment and K is the curvature. If it is assumed 
that the relationship between the curvature and bending moment is linear 
25 (M = k b K), the bending deformation energy can be obtained by the following 

equation. A non-linear case will be described later. 



As illustrated in FIG. 2, prior to buckling before compressive force is 
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(2) 



Where, Kb is the flexural rigidity. 

Since the arc length of the beam is the same as that of the initial 
straight length Li of the beam, the curvature K can be expressed solely in terms 
5 of the distance (| x |) between the two extremities of the beam as follows; 



2 . \ x ij] sin x 
K — — sin c (- — -) y sin c(x) = 



(3) 



The force vector is derived as, 



f t =-^-^k b KL, dK *" =kb K^{^lr^0- (4) 

dx, '<*!*,, II*,, I dK |x,,| 



= ^ 2 (cos^L-sinc(^ Z l))- 1 -^- (5) 

2 2 |x.,| 
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FIG. 5 depicts a relationship between the force vector (compressive force) 
and the distance between the two extremities of the deformation unit. As 
illustrated in FIG. 5, since the immediate buckling model is used in the present 
15 invention, when the compressive force is applied to cloth, the cloth starts to 
buckle at the onset of applying the force (immediately bent without contractive 
deformation). Then, when compressive force exceeding a critical value (Per) is 
applied, a sizable amount of deflection (buckling) occurs. However, in order to 
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simplify the relationship, deformation occurring below the critical value is 
approximated as a straight line as shown in FIG. 6. In other words, 



{ f b , otherwise 

Where, Cb is an arbitrarily selected constant. Although the curvature is 
approximated as described above, there is no significant difference in the results. 
The Jacobian matrix of the force vector thus described is derived as : 

dx t d\x ijX \xl x x m | x] jX x ijX 



In the above equation, — is always positive, but is always 

d I * yl | | | 

negative, such that there is a possibility that the second term of the equation 
10 could cause a singularity to the system matrix. Therefore, the second term is 
dropped. The Jacobian matrix without the negative term, obtained through the 
above-mentioned procedures, is the final result for analyzing the bending 
deformation of the cloth using the immediate buckling model. 

Although in the above-mentioned embodiment, the region at the two 
15 extremities of a predetermined deformation unit of the cloth is modeled as a 
beam having no surface area. The region can be modeled as a plate having a 
surface area. If the region is modeled as a triangular plate as shown in FIG. 3, 
the bending deformation energy can be calculated as follows: 



'• 'J. ,•■"•./.,.;*.,.< (8) 

20 Where Aiand A, respectively define areas of triangles. 



In this case, the force vector can.be obtained as the following equations. 
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Procedures hereafter are the same as those of the beam, except that in 



this case, the geometrical shape of the unit at an equilibrium state is like a 
cylinder. In the above embodiment, although it is assumed that a stable 
structure at the equilibrium state is an arc or cylindrical shape, other shapes 
may be used, since the simulation results are slightly affected by the shape of 
10 the unit at an equilibrium. 

(BENDING DEFORMATION ANALYSIS OF CLOTH BY HYSTERESIS MODEL) 

When the deformation energy is obtained in the aforementioned 
immediate buckling model, the curvature and the bending moments of the 
deformation units are assumed as being linear. Although this assumption can be 
15 well applied to such cloth as silk, other cloth materials such as cotton may not 
be applicable because it has a hysteresis phenomenon. FIG. 7 illustrates a pure 
cotton having hysteresis characteristics in relation with the bending moment and and 

e. 

Conventionally, hysteretic curves obtained by measurement were directly 
20 used to simulate the hysteresis characteristics of cloth. However, the hysteretic 
curves obtained by the measurement can be applied only to a case where the 
curve changes during measurement are identical to curve changes during 
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simulation such that they are not generally appropriate to be applied to dynamic 
simulation. Therefore, it is desirable to make a physical model in order to 
simulate hysteresis characteristics of cloth. 

Meanwhile, generally hysteresis phenomenon is deeply involved with 
5 plasticity. In order to analyze the hysteresis phenomenon, it is widely known to 
use such hysteresis models as shown in FIGs. 8 and 9. FIG. 8 illustrates a slip 
unit having an ideal plasticity to depict a relationship between curvature (X) and 
moment (F). FIG. 9 illustrates a model in which springs and a slip unit are used 
to model a basic hysteresis characteristic. However, the model having the basic 
10 hysteresis characteristics shown in FIG. 9 is rather inappropriate to simulate an 
actual hysteresis characteristic of cloth as shown in FIG. 7. 

Therefore, in order to simulate an intrinsic hysteresis characteristic of 
cloth (different from the basic hysteresis characteristic), it is necessary to use a 
hysteresis model that is appropriate to the cloth. In the present invention, a 
15 new hysteresis model is provided to simulate the actual hysteresis 
characteristics of the cloth as shown in FIG. 7. 

As an embodiment of the present invention, a model is illustrated in FIG. 
10 wherein 4 springs and 3 slip units are connected. The entire responsive curve 
derived from the model is composed of finite straight lines, and the slopes of 
20 each finite straight lines are K0+K1+K2+K3, K0+K2+K3, K0+K3 and K 0 , respectively. 
A more generalized hysteresis model is composed of n + 1 springs and n slip 
unit. 

As illustrated in FIG. 10, when the moment is reduced, the curve follows 
the lines of A4A5, A5A6, AeAy, A?A 8 . The lines of A4A5, A 5 A 6 , A 6 A 7 and A 7 A 8 are 
25 parallel with OAi, A1A2, A2A3 and A3A4, respectively, but the lengths of the lines 
A4A5, AsAe, AeA? and A?As are twice the length of the lines of OAi, A1A2, A2A3 and 
A3A4, respectively. Therefore, asymmetry occurs. 

10 



In the present invention, any hysteresis characteristic of a particular, 
cloth can be expressed by adjusting the numbers of the slip units and springs, 
critical moments of each slip unit and constant values of each spring. As a 
result, the hysteresis characteristic of a particular cloth may be given as a 
5 function of moment M (K) against curvature K. 

In order to apply the above-mentioned new hysteresis model to the 
bending deformation analysis of cloth, it is enough to substitute M in the 
equation 8 for the function M (K). Hereinafter, the above-identified hysteresis 
model is referred to as " Choi-Ko Discrete Hysteresis Model 

10 In FIG. 11, another embodiment of the present invention is illustrated 

which is identical to the finite hysteresis model except that 'n 1 is infinite. The 
curve derived from this model is smooth as shown in the right side of FIG. 11 

This model is used in analyzing the hysteresis characteristics of cloth 
from a point where the plasticity deformation occurs to a point where the 
15 plasticity deformation ends (for example, between Al, A3 and A5, A7 shown in 
FIG. 10) as illustrated in FIG. 11. 

In this case, when the yield limit of the slip unit which reaches the yield 
limit first is given as amin 'and the yield limit of slip unit which reaches the yield 
limit last is given as amax, the yield limits of the infinite numbers of the slip units 
20 are distributed between [amin, amax]. Additionally, it is assumed that the 
(infinitesimal) constants of all springs connected in a series with the slip units 
are the same. 

At this time, according to the descriptions of FIG. 10, when K is given as a 
total sum of the spring constants divided by the length of the portion where 
25 plasticity deformation occurs, as the change rate of the slope at the curve is 



constant, — 
dx 



dX 



K. In this case, if two conditions are given, the second 
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order derivative corresponding, to the curved portion at. the. hysteresis curve can 
be uniquely determined. Considering that the length of the lower portion of the 

d (dF^ 



curve is twice the length of the upper one, it can be modeled as 

dx 



K dXj 



0.5K 



It should be noted that only the absolute value of the second order 
5 derivative corresponding to slip units exceeding the critical values becomes 0.5K, 
and for the slip units not yet exceeding the critical values the absolute value of 
the second order derivative should be K. 

The curve illustrated at the right side of FIG. 11 is an embodiment of a 
hysteresis curve derived from the infinite hysteresis model of the present 
10 invention. 

One of the strengths of the infinite hysteresis model is that, it can reduce 
storage for storing each condition of the slip units and can simulate the 
hysteresis characteristics more realistically in comparison with the finite 
hysteresis model. The infinite hysteresis model thus described is referred to as 
15 "Choi-Ko Continuous Hysteresis Model 

(STRETCH DEFORMATION ANALYSIS OF CLOTH BY SPRING MODEL) 

When stretching force is applied to cloth, the movement of cloth can be 
simulated by spring models installed at two extremities of a predetermined 
deformation unit of cloth. Stretching deformation energy accumulated on the 
20 predetermined deformation unit by the stretching force is obtained to simulate 
the movement of cloth in relation to the stretching. The stretching deformation 
energy can be derived from the following equation. 
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E = 



(10) 



x ij2 H L 2 



where x ij2 is the distance between the two extremities of the 
predetermined deformation unit, L2 is the distance where the stretching force is 
not applied, and Ks is a spring constant. In case of \x lj2 \<L 2 , the deformation 

energy is 0, such that it should be noted that the above model can be applied to 
cases where stretching forces are applied among the extremities. In the above 
equation, the force acting on a certain position between the two extremities of 
the predetermined deformation unit can be expressed as follows in Equation 11: 



0, I x.. H L 2 



(11) 



Where, the force vector of the Jacobian matrix is given by Equation 12. 
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The first term of the above Jacobian matrix is derived from modeling the 
interaction as a linear spring and it indicates that the stiffness along the 
interacting direction by the stretching is constant. The second term of the 
Jacobian matrix indicates that the stiffness along the perpendicular direction of 
the interaction is proportionate to (1-—^—) direction. Since the cloth is shaped 

I *V2 I 

in a planar form, the perpendicular direction of the interaction corresponds to a 
direction perpendicular to the plane. 
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Because the spring model is applied only to a case where stretching force 
is applied, |x | is larger than L 2 , and is positive but smaller than 1, 

1**2 I 

thereby rendering the system stable. If the spring model is used in compression 
case, as \x 2 \ approaches to 0, (1--^—) diverges to infmiteness. Therefore, 

5 the entire system matrix (I -a — + p — ) may generate a singularity problem 

dv dx 

regardless of the simulation step size. Therefore, the spring model is not used 
for compression in the present invention. 

By way of reference, in the prior art [ Volino and Magnenat-Thalmann 
2000 ], the same spring models are used in both stretching and compression, 
and the second term of the Jacobian matrix in Equation 13 is dropped in both 
the stretching and compression. Furthermore, in order to avoid eigenvalues in 
the orthogonal direction of the interaction from becoming null due to the drop of 
the second term, a fictitious damping term is added. Although the the adde 
damping term can stabilize the system matrix, it may unnecessarily generate 
high resistance to the movement of the cloth. 

As for compression in the present invention, an immediate buckling 
model is used which causes an instant bending deformation without contraction, 
while as for stretching, a spring model is used, such that stability of the system 
matrix can be obtained without adding an imaginary damping term. Therefore, 
20 the present invention can simulate the movement of cloth more realistically. 

As apparent from the foregoing, there is an advantage in the immediate 
buckling model, hysteresis model, and simulation method based on the invented 
models and computer-readable media that stores the program which executes 
the invented simulation method thus described according to the present 
25 invention in that the immediate buckling model is used to simulate deformation 
caused by compressive forces, thereby solving the structural instability of the 
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buckling phenomenon. Additionally, the simulation speed is increased. There 
is another advantage in that " Choi-Ko Discrete Hysteresis Model "and " Choi-Ko 
Continuous Hysteresis Model " are used to simulate the movement of cloth, 
thereby simulating the hysteresis characteristics more realistically. 

While the invention has been described in terms of its preferred 
embodiments, the invention is illustrative and not restrictive thereto, and may be 
embodied in several forms without departing from the spirit or essential 
characteristics thereof, and it should be apparent that those skilled in the art will 
recognize with considerable modification within the spirit and scope of the 
appended claims. 



15 



